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Abstract 

We present a modular cosmology scenario where the difficulties encountered in conven- 
tional modular cosmology are solved in a self-consistent manner, with definite predictions to 
be tested by observation. Notably, the difficulty of the dilaton finding its way to a precarious 
weak coupling minimum is made irrelevant by having eternal modular inflation at the vacuum 
supersymmetry breaking scale after the dilaton is stabilised. Neither this eternal inflation 
nor the subsequent non-slow-roll modular inflation destabilise the dilaton from its precarious 
minimum due to the low energy scale of the inflation and consequent small back reaction on 
the dilaton potential. The observed flat CMB spectrum is obtained from fluctuations in the 
angular component of a modulus near a symmetric point, which are hugely magnified by the 
roll down of the modulus to Planckian values, allowing them to dominate the final curvature 
perturbation. We also give precise calculations of the spectral index and its running. 

PACS: 98.80.Cq 



1 Introduction 

Early Universe cosmology offers one of the few opportunities to test the framework of particle theory 
phenomenology at high energy scales, in particular string/M-theory phenomenology. Indeed, moduli 
are generic in string/M-theory, and much attention has been drawn to cosmology where moduli 
play important roles in the early Universe pp. This has been referred to as 'modular cosmology'. 
Unfortunately, despite many attempts to find viable cosmology scenarios associated with moduli, 
they turned out to give us cosmological obstacles |21 El EI rather than attractive features. 

Our main goal in this paper is to present a successful modular cosmology scenario where those 
difficulties due to moduli can be solved in a simple and self-consistent manner. Furthermore, this 
scenario gives definite predictions to be tested by forthcoming observations. 
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The paper is structured as follows. 

Section 12.11 gives our motivation for the investigation of modular cosmology and discusses the 
cosmological difficulties involving moduli. Section [2.21 outlines how those difficulties are overcome 
in our proposed scenario. 

Section IH.ll discusses eternal inflation at the vacuum supersymmetry breaking scale, which 
is crucial in solving the dilaton stabilization problem. Section 13.21 is devoted to the non-slow-roll 
modular inflation which naturally follows the eternal inflation, and whose low energy scale is crucial 
not to disturb the dilaton potential. 

Section IH gives a detailed analysis of the density perturbations produced in our modular cos- 
mology, starting in Section 14.11 with a brief review of the 5N formalism jH] necessary to analyse 
the contributions of the fluctuations in the multiple light degrees of freedom to the final curvature 
perturbation. 

Section El briefly discusses thermal inflation's solution of the Polonyi/moduli problem. 
Section ini gives our conclusions. 

Facing the fact that our knowledge of string/M-theory is still too limited to fully describe 
the moduli dynamics in early Universe, we keep our discussion on string/M-theory aspects fairly 
general with emphasis on cosmological consequences to be compared with data. We hope our 
scenario can open up new possibilities in model building of modular cosmology which can be tested 
by observations. 

2 Motivation 

2.1 Difficulties in Modular Cosmology 

Moduli are generic predictions of string/M theory, but pose major cosmological problems. 

The stabilization of the dilaton is a serious problem in both particle physics and cosmology. In a 
typical setup, one simply cannot find a minimum in the weak coupling regime [Hj. More complicated 
models, such as the racetrack scheme [Zj, can realize a minimum at weak coupling with a moderate 
amount of fine-tuning (which is easily anthropically justified if there is no better solution). However, 
even if there exists such a minimum, it has only a shallow barrier separating it from the runaway 
part of the dilaton's potential, and a steeply rising potential on the other side. Thus any typical 
cosmological history would lead the dilaton to over-shoot the desired weak coupling minimum and 
roll off along the runaway potential towards zero coupling |3] . 

Another serious cosmological problem is the cosmological moduli problem. Many moduli are 
expected to have Planckian vacuum expectation values relative to any symmetric points where they 
couple to other fields, giving them gravitationally suppressed couplings to other light fields. Shifting 
of their minima due to supersymmetry breaking in the early universe causes them to dominate the 
energy density of Universe, and they survive until or beyond nucleosynthesis destroying its successful 
predictions 0. 

Thus, the ubiquitous moduli predicted by string/M-theory have given us a cosmological headache 
rather than helping us build a desirable cosmology scenario. It is therefore worth while to seek 
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realistic modular cosmology scenarios where those serious cosmological problems are solved without 
fine-tuning in a self-consistent manner, and we address such a possibility in this paper. 

2.2 Successful Modular Cosmology 

Let us outline the key ingredients in our successful modular cosmology scenario. 

One of the main features in our scenario is eternal modular inflation at the vacuum supersymme- 
try breaking scale. This is important because the inflnite expansion of eternal inflation compensates 
for any flnite unnaturalness in initial conditions. In particular, it can compensate for the extreme 
fine tuning of initial conditions required for the dilaton to find its way into a precarious weak cou- 
pling minimum. The importance of the low energy scale is that it ensures that the dilaton does not 
get disturbed from its minimum after the eternal infiation ends locally. If one had eternal inflation 
at a higher energy scale than the scale of vacuum supersymmetry breaking, the backreaction on the 
dilaton's potential would dominate over the dilaton's vacuum potential, and one could never expect 
the dilaton to flnd its way into its desired weak coupling minimum after the eternal inflation ends. 
Moreover, if eternal inflation occurs at a higher energy scale than the vacuum moduli potential, the 
moduli will roll down to possibly several vacua including run-away vacua, also leading to domain 
wall problems. 

Of course, we do not mean that all of eternal inflation takes place at this low energy scale, just 
that the last eternal inflation (from our local point of view) occurred at that low energy scale. The 
full picture of eternal inflation would be of eternal inflation at the various discrete points in fleld 
space, from the Planck density down, corresponding to suitable false vacua or maxima with rare, but 
finitely rare, transitions between these points, and of course continually generating non-eternally 
inflating regions like our own [8| El ^] • 

Since the last eternal inflation occurred at the vacuum supersymmetry breaking scale, the in- 
flation that produced the observed approximately scale-invariant density perturbations must have 
occurred at that scale or below. Such inflation emerges naturally as the modulus rolls away from its 
eternally inflating maximum. The scale-invariance is not automatic but can arise naturally when 
proper account of the contributions of the fluctuations in all the relevant light degrees of freedom 
in the generically complicated moduli potential is taken into account. In particular, fluctuations in 
angular degrees of freedom are hugely magnifled by the long roll down of the modulus to Planckian 
values, allowing them to dominate the flnal curvature/density perturbation. 

Finally, the overabundant moduli are diluted by thermal inflation which is natural consequence 
of thermal effects on unstable supersymmetric flat directions [llj. 

We now look into the details of our scenario in the following sections. 
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3 Eternal and Post-Eternal Modular Inflation 



3.1 Eternal Modular Inflation 

The modulus potential has the general form 

V{^) = M^r{^lM^{) , (1) 



where Ms is the scale of supersymmetry breaking and Mpi = 1/ vSttG. $ is the scalar component 
of a modulus chiral superfield and JF is a function with dimensionless coefficients of order unity. 
For most of the rest of the paper, we set Mpi = 1. In our vacuum, for the case of gravity mediated 
SUSY breaking, Ms ~ 10^° — 10^^ GeV, in order to get supersymmetry breaking masses rris ~ ~ 
10^ — 10^ GeV. In our regime of interest, these will be the relevant scales. ^ 
Near a maximum, this potential will have the form 

V(^) = Vo- \ml4>' + . . . , (2) 

where Vq ~ and ~ M^^, and we have simplified to a single real direction G $. The equation 
of motion is 

+ 3iJ0 - ml<i) = (3) 
which, in the approximation 3H^ ^ Vq, has growing solution 

oc e"^* oc a" , (4) 

where 

^ = ^\ + ^ -1 I • (5) 




Classically, one will get eternal infiation at the maximum as the field never rolls off. Quan- 
tum mechanically, one will have a distribution of field values which will scale like a". Thus, the 
probability of having a field value |(/)| < e for some small e will go like 

P(|0|<e)oc-^cxa-^ (6) 

while the physical volume of space scales as a^. Thus if a < 3 or 

ml < QVo , (7) 

one will have eternal infiation at the maximum [Tn| I13j. 

As we discussed in Section 12. 2| eternal infiation at this low energy scale plays a crucial role for 
stabilization of the dilaton in our scenario. 



^We expect our scenario to work equally well in the case of gauge mediated supersymmetry breaking, but leave 
proper investigation to |12| . 
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3.2 Post-Eternal Modular Inflation 



Once the modulus locally rolls away from the maximum, the eternal inflation will end locally, but 
one will still get some more inflation as it rolls towards its minimum at Planckian values. This 
inflation is important because fluctuations produced during the eternal inflation are large leading 
to an inhomogeneous universe. The spectrum of perturbations produced by the modulus as it rolls 
away from the maximum is 



r (I) 27ra0^ \a^H, 

where subscript -k denotes an appropriate evaluation point ^. For our values of ~ 10^^^ and 
a ~ 1, this drops below the COBE normalisation for 

> 0COBE ~ 10-'' . (9) 

Thus to make the observable universe sufficiently homogeneous we require that observable scales 
leave the horizon some time after cf) = 0cobe- The observationally relevant number of e-folds is the 
number of e-folds from this point until inflation ends at </> ~ 1, which is 

N = Un(-^)^^-^. (10) 
" V0COBE/ a 

Recalling that, for inflation at l/^/^ ~ 10^° - 10" GeV and assuming radiation domination after the 
inflation, the observable universe leaves the horizon around 45 e-folds before the end of inflation 
|15j . and taking account of an additional ~ 10 e-folds due to thermal inflation, we see that for a ~ 1 
we have enough e-folds of inflation to make the observable universe sufficiently homogeneous. 



4 Density Perturbations 

Beyond making the universe homogeneous, we also need to produce the observed approximately 
scale-invariant spectrum of density perturbations. The previous simplified model fails to do this for 
the usual reason: the supersymmetry breaking effect of the inflationary potential energy density 
induces a mass squared m| ~ Vq and hence we have a ~ 1 and a spectrum far from scale invariant 

uniiizi. 

However, to properly calculate the spectrum of density perturbations produced during modular 
inflation, we need to take proper account of the contributions of fluctuations in all the relevant light 
degrees of freedom to the flnal curvature/density perturbation, and not just use an over-simplifled 
toy model. We start with a brief review of the formalism for this jH] . 

^One would usually choose some time around horizon crossing as the appropriate evaluation point, but this is not 
necessary here. In fact, the evaluation point can be any time during inflation when the scaling of 0, (f> cc a", and the 
constancy of H ensure that Eq.(8) is independent of the time chosen to evaluate it. 
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4.1 SN Formalism 



The 6N formalism jB] relates the final curvature perturbation to the perturbation in the number of 
e-folds of expansion from an initial spatially flat hypersurface to a final comoving, or equivalently 
constant energy density, hypersurface 

n^ = 5N . (11) 

For illustration, let us consider single field dynamics first. In this case, the field trajectory is 
unique and therefore any perturbation just corresponds to a shift along the trajectory, which can be 
compensated by a local time shift, and hence leads to a curvature perturbation depending only on 
local quantities. If there exist multiple fields, however, there is an infinite family of field trajectories 
in the multiple-field space and each of the trajectories can follow different equations of state. Thus 
perturbations which shift from one trajectory to another can change the equation of state history 
and hence the number of e-folds of expansion 

N= [Hdt = - [h^= (12) 
J J p J 3{l + w) ^ ' 

in a very non-local way. Here w = p/ p. The final curvature perturbation is achieved when the mul- 
tiple field dynamics disappears and the trajectories coalesce, leading to constant TZc on superhorizon 
scales. 

To apply this formalism, we simply write ^ 

ON 

n, = 6N=—6(j),, (13) 

where 0j represents each dynamical degree of freedom. Note that, in the above expression, dN/dipi 
is a transfer function depending on background fields while Sipi is the perturbation on the initial 
flat hypersurface, and they can be calculated separately. 

In the following sections we apply the 5N formalism to calculate the final curvature perturbation 
in three concrete toy models properly taking into account all relevant degrees of freedom. ^ 

4.2 Simple Single Modulus 

As an introduction, let us discuss the simplest case of a single complex modulus 

$ = -L^e^^ . (14) 



Eq. dni) gives 



dN dN 

n^^m^ + (i5) 



In general one has SN = §^ S(f>, + 1^ S^i but usually one can pick out a growing mode for S(f>i in which case 

S(l)i and S(j)i are no longer independent. 

''in this paper we mostly apply the 5N formalism to dynamics towards the end of inflation. However, it also 
encompasses cases, such as curvaton scenarios jl8j . where the multiple degrees of freedom survive long after inflation. 
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We assume that $ = is a symmetric point so that the potential near the maximum has the form 



V = VQ-ml\^^ + ... . (16) 

However, now that we are considering multiple real fields, just knowing the local potential is not 
good enough. The full modulus potential will be completely model dependent but a simple toy 
model will be sufficient for our purposes. We take 

V = Vq- mj|$|2 + \ml {^^ + h.c.) + mjl^l^ (17) 

with Vq = 2m^/3 for vanishing vacuum energy. 

The contribution to the final curvature perturbation from the radial component is 

^^i^^H^^" (18) 

and was given precisely in Eq. (jH)). The contribution to the final curvature perturbation from the 
angular component is 

To evaluate this we need to determine dN/d9. 

The analytical estimate of ON/ 89 in our scenario is non-trivial because the approximate U(l) 
symmetry at small is heavily broken at large 0. We hence give numerical analysis to determine 
ON /do. As our initial conditions, we take a small value of 0, where the potential is still angularly 
symmetric so that the exact value of does not matter, and various values of 9. Some field 
trajectories are shown in Figure ^ 

We can see that the initial angle 9i greatly affects the form of the trajectory for large field values 
and will lead to different equation of state histories, and hence different amounts of expansion. 
A'"(^) is shown in Figure 2, where we integrated the equation of motion starting from a fixed small 
initial value of and continuing until a fixed final energy density where the equations of state had 
converged ^. 

We see that 

dN 

-QQ - Oil) (20) 

for a wide range of 9i, though the initial values 9i ~ 2n7r/3, which lead to large spikes in N, could 
be a postiori preferred due to the greater spatial volume produced [T^ . 

Thus, from Eqs. (fTH|) and (fT^ . we see that the angular fiuctuations give a comparable, or possibly 
greater, contribution to the final curvature perturbation. This is due to the huge magnification of 
the angular fiuctuations as rolls from small values to Planckian values, comparable to the usual 
huge magnification of radial fiuctuations near a maximum due to the factor dN/dcj) = —H/cf). 

^Here we ignore complications present in more complete models where one has to take into account degrees of 
freedom suppressed here. In particular, one might expect the different trajectories to lead to different fractions 
of radiation or other things produced at the end of the modular inflation. We will discuss this in more detail in 
Section 
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Figure 1: Trajectory of $ = ('I'r + i$i)/-\/2 for initial angles vr/Q and 7r/3. 

However, as dN/dO does not depend on the initial value of (p (assuming that it is small), dN/d6 
will be independent of k. Thus the contribution of the angular fluctuations, Eq. (fT^ . has the same 
steep spectrum as that of the radial fluctuations and is not consistent with observations. 

In the next sections we introduce models which avoid this problem. 

4.3 Single Modulus with Quantum Corrections 

Here we consider the next simplest model of modular inflation given in Ref. ^H]- We first briefly 
describe the model and then discuss the perturbations produced by it. 

The previous model considered a modulus rolling away from a maximum at a point of enhanced 
symmetry. At such a point of enhanced symmetry, one would in general expect the modulus to have 
couplings to other light fields. These couplings will cause the mass to renormalise as a function of 
(j). This could turn the mass squared of the modulus positive for small 0, shifting the maximum of 
the potential out to some finite value = 0^,. 

Explicitly, following Ref. ^Hl with the notation change e ^ (3 and implicitly keeping the full 
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Figure 2: Number of e-folds as a function of the initial angle 9i. 
complex modulus field $, 



V = V^ 

and we have a maximum aX (h = (h.^ where 



l-^/(/?ln</.) + 



(21) 



/* + \^fl = . (22) 

We resolve the ambiguity in the definitions of / and P by fixing = /(—I) so that 0^, = e~^^^. 
Note that /* < and typically /3 -C 1. We assume H ^ a/Vo/3 and 0^0* when observable scales 
leave the horizon. 

One will again have eternal inflation at the maximum, which now corresponds to the ring cf) = cf)^, 
and the eternal inflation ends locally when locally rolls off the maximum towards large values. 
The distribution of 9 values during the eternal inflation could be biased by the small symmetry 
breaking effects of the terms, which should also be included in this model, because the difference 
in expansion rate around the ring leads to differences in the spatial volumes generated. However, 
for small 0*, this biasing will be negligible as the differential expansion rate is much smaller than 
the diffusion rate around the ring due to quantum fluctuations. 

This model has the interesting property that in the neighborhood of the maximum the potential 
is flattened in the radial direction by a factor (3. This helps to flatten the contribution of the radial 
fluctuations to the flnal curvature perturbation, but not by as much as one would like jTH] . However, 
we will be more interested in the effect of the shift of the maximum out to flnite = 0=,,. This means 
that as rolls off the maximum, while — 0* ^ 0*, the radius will hardly change, allowing the 
angular fluctuations to give a flat contribution to the flnal curvature perturbation. 

Explicitly, the contribution of the radial fluctuations to flnal curvature perturbation is 
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and the contribution of the angular fluctuations is 

dN dNH 



09 89 



(24) 



Therefore, for dN/d9 ~ 0{1), the contribution of the angular fluctuations is suppressed relative to 
that of the radial fluctuations by the factor (0 — 0*)/0. This is easily understood as for large 0, 
i.e. on small scales, the contributions are comparable, but as one approaches 0*, i.e. going to larger 
scales, the angular contribution to the spectrum flattens out while the radial contribution continues 
rising and so dominates. 

However, if ON/ 89 3> 1, then the contribution of the angular fluctuations can dominate the 
final curvature perturbation giving the desired flat spectrum. ^ This is realized, as indicated from 
Figure 121 if the angle is close to the one of angular maxima of the potential at ^ = 2n7r/3. As 
mentioned before, this can be probable if we consider that regions corresponding to these initial 
angles will undergo greater expansion and hence occupy a larger volume at late times. However, we 
do not want to be drawn too close to these special angles as dN/d9 will either diverge or be zero 
there. This requires further investigation |l2j. 

A more precise formula for the spectrum of 69 is given by Eqs. and in the Appendix. 
The spectral index has the form 

n-l = -Ak^f' , (25) 

and the running 

Therefore, if the contribution of the angular fluctuations dominates, one can expect a flat spectrum 
on larger scales with the spectrum dipping away with substantial running on smaller scales. It is 
very plausible that observable scales correspond to (p close to, but not very close to, 0*, so that this 
running would be observable |2Uj . 



4.4 A Second Modulus 

String theory does not predict a single modulus, but rather many moduli. Here, in addition to 
we consider the case that there is a second modulus \E' which drives eternal modular inflation and 
also some post-eternal modular inflation, before $ rolls down driving its own modular inflation. 
We take \&'s potential to have the form 

V = Vo- \m\^^ + . . . , (27) 

where we consider just a single real direction t/^ G and, as always, ~ Vq and everything is at 
the vacuum supersymmetry breaking scale. As before, ■?/' rolls down as 

oc a" , (28) 
rough estimate suggests we need dN/dO > 10^. See Appendix IA.2I 
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where 

and the contribution of fluctuations in ijj to the final curvature perturbation will be 

dip ip ip 

During the modular inflation driven by ip, $'s potential will have the form 

V = ViiP) + ^ml{<P,iP)^' + ... . (31) 

The mass squared rn^lcpjip) is renormalised as a function of as before. However, it will also have 
a dependence on ip due to the supersymmetry breaking effect of ip^s potential. Schematically, it 
looks like 

ml{(j),ip) ml^^ + mj;^ + /5m^^^ In , (32) 

where m'^^^ ~ mf^^ ~ m^g^ ~ and /3 -C 1. We assume the indicated signs and that mf^^ > m1^^ 
so that > at large during -j/^'s modular inflation. We further assume that the renormalisation 
group running eventually turns negative so that has a minimum at some exponentially small 
value 0*. This is just the inverse of the running we had in the previous model which produced a 
maximum, and leads to a quite interesting scenario as follows. 

During -^'s modular inflation, has a mass of order H which will suppress its fluctuations 
while ^'s angular component stays massless and so will have unsuppressed fluctuations. Once t/^'s 
modular inflation ends, m^^j will disappear, or at least change greatly, releasing to roll down to 
its vacuum at Planckian values as in our previous models, driving some more modular inflation 
as it does so. The angular fluctuations in $ are massively amplified by this roll down and give a 
contribution to the final curvature perturbation 

dN ON H 

assuming that observable scales leave the horizon during ip^s modular inflation. Since ON/ 89 is 
independent of /c, and 0* is a constant, we get a flat contribution to the final curvature perturbation. 
Furthermore, it will dominate over -^'s contribution as long as -0 ^ 0* when observable scales 
leave the horizon. Finally, for H ~ 10~^^ and ON/ 89 ~ 0(1), our model matches the COBE 
normalisation with a very reasonable 0* ~ 10~^^, which is a natural consequence of the logarithmic 
renormalisation group running in our weak coupling regime. 

A more precise formula for the spectrum of 59 is given by Eqs. and (|Un|) in the Appendix. 
The spectral index has the form 

n - 1 = -Afc^a ^ ^34) 

^Note that the flattening mechanism also apphes in this case, so the fluctuations in (j) may not be completely 
suppressed. We will leave further investigation of this point to another paper 
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and the running 



'^'^ -2aAA;2- . (35) 



d\nk 

So, as in our previous model, one can expect a flat spectrum on larger scales with the deviation 
from a flat spectrum on smaller scales. In this model one expects a ~ 1 so the running will typically 
be larger than in the previous model. Due to the second modular inflation driven by 0, it is very 
plausible that observable scales left the horizon near the end of the first modular inflation driven 
hy ip, so that this running may be observable |2Uj . 

As mentioned in the Appendix, the above formula for the spectral index applies only for a < 1. 
The case a > 1 is more involved because the deviation from scale invariance is dominated by 
non-trivial dynamics at the end of inflation. Thus, in our model, even though the spectrum will be 
rather flat on large scales, we could possibly expect deviation from a flat spectrum for smaller scales 
caused by late-time inflation dynamics affecting super-horizon modes ^]. We leave the detailed 
analysis of this perturbation calculation on super-horizon scales during the late epoch of modular 
inflation for our forthcoming paper 



5 Thermal Inflation 

After our modular inflation, the moduli roll down and oscillate about the minimum of their potential. 
If this minimum corresponds to a point of enhanced symmetry, where the modulus has unsuppressed 
couplings to other light fields, then the modulus will decay rapidly. If the minimum is a Planckian 
distance away from any such point, then the modulus will have only gravitational strength couplings 
to other light fields and will decay too late for the successful predictions of Big Bang nucleosynthesis 
to survive [2]. Even if our modulus decays rapidly, one expects other moduli to be displaced during 
the modular inflation and to roll down and oscillate when it ends. Again, these moduli may either 
decay rapidly, or contribute to the Polonyi problem. Thus, after the modular inflation, one expects 
to have initially perhaps comparable fractions of radiation and dangerous long-lived moduli. 

Thermal inflation is considered to be one of the most promising mechanisms to get rid of the 
unwanted moduli. We keep the following discussion of thermal inflation fairly brief and refer the 
reader to Ref. JT] for more details. 

Thermal inflation is driven by a flaton a with potential of the form 

V = Vo+ [g^T^ - m^) \a\^ + {AW + h.c.) + |A| Vl'""' , (36) 

with m ~ 10^ — 10^ GeV and n = 4 or 5. The flaton's vacuum expectation value is M ~ ■mi/C"-^) 
and Vo ~ m^M^ -C so we do not expect eternal or modular inflation. When T > Tc = m/g, 
the flaton is trapped at the origin by the finite temperature potential, and inflation begins when 
Vo dominates over the moduli density. After about 5 to 10 e-folds of inflation, which dilutes the 
moduli, the temperature drops below and the flaton rolls rapidly towards its minimum ending 
the inflation. The Universe is then dominated by the oscillating flaton, which eventually decays 
to ordinary fields through non-renormalisable operators suppressed by of order M. After this 
'reheating' at temperatures of order GeV, standard Big Bang cosmology follows. 
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We point out here that fluctuations in the angular part of the modulus produced during the 
modular inflation, which lead to perturbed trajectories at the end of the modular inflation, could 
lead to perturbations in the fraction of radiation produced at the end of the modular inflation. This 
could then lead to perturbations in the amount of thermal inflation, as the thermal inflation begins 
when the moduli density drops below Vq but ends when the temperature drops below Tc, giving an 
extra contribution to dN/dO. However, the perturbation in the fraction of radiation produced at 
the end of the modular inflation is likely to be highly model dependent and difficult to evaluate. 
This will be further investigated in Ref. [T2|, though we do not expect it to change our conclusions 
greatly. 

We also note in Appendix IA.3I that thermal fluctuations can produce perturbations during 
thermal inflation, but these are negligible in our case. 

Thus thermal inflation, which is a natural consequence of thermal effects in the presence of 
supersymmetric flat directions, resolves our remaining difficulty with modular cosmology. 

6 Conclusions 

We presented the first realistic modular cosmology scenario in that dilaton stabilization, inflation 
producing a scale invariant spectrum of density perturbations, and Polonyi/moduli problems are 
resolved in a simple and self-consistent manner. 

One of the notable features is eternal inflation at the vacuum supersymmetry breaking scale. 
This causes the fine-tuned initial conditions required for the dilaton to find its way to a weak 
coupling minimum to become a postiori as natural as any others because of the subsequent eternal 
inflation. 

Our model also predicts the observed scale-invariant cosmic microwave background spectrum 
without any fine-tuning, and even allows for the possibility of a dipping spectrum with significant 
running on smaller scales [20]; which alone could deserve special attention from the viewpoint of 
inflation model building. Indeed, non-slow-rolling inflaton fields are generic in supergravity theory, 
but usually avoided by fine-tuning because of their non-scale-invariant power spectrum. This fine- 
tuning is not necessary in our scenario because the density perturbation is dominantly produced 
from the angular degree of freedom of a modulus whose usually subdominant fluctuations are 
hugely amplified by the long roll out of the modulus. We gave explicit calculations of the spectrum, 
spectral index and its running for modular inflation for concrete toy models. We also pointed out 
the possible deviation from scale invariance at small scales for some values of the parameters, and 
we leave thorough treatment of this issue for our forthcoming paper [T2] . 

Finally, thermal inflation solves the Polonyi/moduli problem. 

Due to the generic existence of ubiquitous moduli in string/M-theory, an early universe history 
where moduli play crucial roles would be a rather natural consequence. We kept our discussions 
fairly general so that it can be apphed to a wide class of models, and hope our proposed scenario 
leaves room for the further investigation of successful modular cosmology. 
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A Appendices 

A.l Precise Calculation of 56 

In this appendix we calculate the spectra for the models presented in Sections 14. HI and 14. 4| starting 
with a general formula applicable to both. 

To calculate 59, on superhorizon scales but while the angular potential is still flat, we need to 
solve the perturbed equation of motion for 69{k,t) 



dt 



3 J,2 



d69 



ak'^(t?59 = . 



Defining 



Lf = a(l)59 

and using the conformal time drj = dt/a, we obtain 



drf 



+ 



1 d'^{a(, 
acj) drf 



= 



where 



1 d^i^ac, 
acf) di]'^ 



1 dH 



+ 



1 d^ 



2iy2 dt 2H(t) dt 2i720 df^ 



Introducing y = y2k(f and x = —krj, we can write this as 



d^ + \'--)y = -y 



x^ 



where 



xacf) 



^ d^{xa(j)) ^d{xa(f)) 
{dlnxy dlnx 



(37) 
(38) 

(39) 
(40) 



(41) 



(42) 
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For 

with ^ <^ 1, this can be solved to give 

2 

Pseik) 



9 = i 



X 



(3 - v)v 



{2x^) ^ cos 



/TH/X V{2 + u) _ ^ 
V 2 / l-v 



(43) 



(44) 



where subscript -k denotes an appropriate evaluation point. The divergence at u = —2 corresponds 
to the deviation of the spectrum from scale invariance being dominated by late time dynamics |21j . 
The spectral index is 



nseik) - 1 = - cos(^— j — 



7Tu\ 2T{2 + u) ^ 



(i-/^)(3-z/) i2x,y ■ 

We apply the above formula to the specific models we discussed in the paper. 

A. 1.1 Single Modulus with Quantum Corrections 

Here we derive the explicit form for Eq. of Section 14.31 
For H ~ a/Vo/3 , we have 

- A 

and 

3 d(j) 1 rfV 

9 



Eq. (24) of Ref. ^ gives 



(pH dt (pm df^ 



-/3fUN+C) 



where N is the number of e- folds until inflation ends at ~ 1, and C ~ 1//?. Therefore 

1 dc 



fie 



and 



Hence 



1 rfV 



l^-f}fl{N+C) 



fie 



(pH dt 



9 = {^+ m /:e-/^/*(^+^) = (3 + (3f:) pf: in ( ^ ) ( ^ 



leading to 



21x6 



1 - 2 



2/^/:cosMl(^^r^ 



In 



(45) 



(46) 
(47) 

(48) 

(49) 
(50) 
(51) 

(52) 
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and 



ns0{k) - 1 = -2^-^* cos 



nPf:\ T{2-Pfi 



2 J {1 + Pfi 

Note the sign change from one's naive expectation for f3fl > 1 [T^ . 



(53) 



A. 1.2 A Second Modulus 

Now we derive the exphcit form for Eq. (jMj) of Section 14.41 

Here we assume (j) = (j)^ leaving consideration of the effects of any time dependence of (p to 
another paper ijj scales as 



with 



a 



3 / 



3Vr 







Therefore 



and 



Hence 



leading to 



X 



aH 



1 + 



2(1 -2a) 



1 /3 + 2a 



2 ; 2 

2 V 1 - 2a y ^'^ \ X, 



X 



-2a 



a 



2„ , ,r(2-2a) / k 
22"cos vra ^ ' ' 



1 - 4a2 V a^H- 



2a 



and 



nsg{k) — 1 = — 2 ° cos(7ra) 



r(2-2a) 
1 -4a2 



2„/.2 



a ip. 



2a 



(54) 
(55) 

(56) 
(57) 
(58) 



(59) 



(60) 



A. 2 Estimation of the Required Value of dN/d9 for the Single Modulus 
with Quantum Corrections Model 

In the single modulus with quantum corrections model of Section 14 ■3[ the radial fluctuations have 
only a slightly flattened spectrum. Therefore we require 



(61) 
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to avoid the radial fluctuations giving the full spectrum too big a deviation from scale invariance. 
Now 

and using the results of App endix I A . 1 . 1 1 we have 

' ~/3/:inf^V^:S10-. (63) 



Thus we get the estimate 



We also note that matching the COBE normalization requires 

dN H 



86 0, 

and so for H ~ 10~^^ require 



10"* , (65) 



^ - l^''ct>. (66) 
which gives us an upper limit on dN/dO as we expect 0* to be well below the Planck scale. 

A. 3 Perturbations from Thermal Fluctuations during Thermal Infla- 
tion 

In this appendix we note that there can arise another source of perturbations caused by thermal 
fluctuations during thermal inflation. These thermal fluctuations can be estimated by considering 
the thermal fluctuation 5T ~ T in regions of the thermal bath of size ~ 1/T. A Hubble volume 
of radius 1/ outside of which the thermal fluctuation freezes in, encloses T^/H^ of such regions 
giving a thermal fluctuation on the Hubble scale of 6T ~ T/ y^T^/H^. Thus we have 

Now T oc a and so this has a spectral index n = 4, with its maximum value ~ ^j^g 
end of thermal inflation. This is completely negligible for thermal inflation with an intermediate 
scale M as we are considering, but could be signiflcant for thermal inflation driven by a flaton with 
M > Mgut- 
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